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Abstract 

Direct CP violations in B14 decays (B^ —* 7r + 7r _ / ± t'/) are investigated within the 
Standard Model (SM) and also in its extensions. In the decay processes, we include 
various excited states as intermediate states decaying to the final hadrons, ir + + 7r~. 
The CP violation within the SM is induced by the interferences between intermediate 
resonances with different quark flavors. As extensions of the SM, we consider CP 
violations implemented through complex scalar- fermion couplings in the multi-Higgs 
doublet model and the scalar-leptoquark models. We calculate the CP-odd rate 
asymmetry and the optimal asymmetry. We find that the optimal asymmetry can be 
measured at la level with about 10 9 i?-meson pairs in the SM case and 10 3 -10 7 pairs 
in the extended model case, for maximally- allowed values of CP-odd parameters in 
each case. 
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1 Introduction 



Semileptonic 4-body decays of 5-mesons with emission of a single pion have been studied 
in detail by many authors JI], ^ Recently we investigated the possibility of probing 
direct CP violation in the decay B ± — > Dtt^u || in extensions of the Standard Model 
(SM), where we extended the weak charged current by including a scalar-exchange inter- 
action with a complex coupling, and considered as specific models the multi-Higgs doublet 
(MHD) model and the scalar-leptoquark (SLQ) models. In the present work, we inves- 
tigate the same possibility in the decay of B ± — > {ti + 'k^)1 ± v. In this case we find there 
may be direct CP violation even within the SM. 

As is well known, in order to observe direct CP violation effects, there should ex- 
ist interferences not only through weak CP-violating phases but also with different CP- 
conserving strong phases. In the decay of — > titt^u, we consider it as a two-stage 
process: B — > Mi — > ttit)Ii>, where M« stands for an intermediate state which is decay- 
ing to 7r + + ix~ . In this picture the CP-conserving phases may come from the absorptive 
parts of the intermediate resonances. Here we try to include as many as possible inter- 
mediate states decaying to 7r + + 7r~, so that they could represent a pseudo complete set 
of the relevant decay. The candidates in b — > u transition are p, fo and f 2 mesons, which 
decay dominantly to 2tt mode (See Table. 1). Furthermore, we find that even in b — > c 
transition a D° meson can decay to tc + + tt~ , although its branching fraction is very small 
compared to those of uu states. However, we can find that the contribution through an 
intermediate D meson is not negligible because of CKM favored nature of b — > c transition 
compared to the b —>■ u one of uu states, p, fo and f2 mesons. If we include D meson as 
an intermediate state as well as the uu states, direct CP violation may arise even within 
the SM through their relative weak phases of the different CKM matrix elements (V c b and 
V u b). Therefore, we first consider CP violation within the SM by including p, f , f 2 and 
D mesons^ as intermediate states decaying to 7r + 7r~. Next we also consider CP violations 
in extensions of the SM, in which we use a cutoff to the final state irn invariant mass so 
that the effects of D meson cannot enter, thus ensuring that the result is solely from new 

° Here we are including fully known resonances only, and neglecting possible non-resonant 2tt decays. 
A significant experimental enhancement can be made, if wc use the reduced kinematic region around 1.4 
GeV < y/(p„ + p„) 2 < 1.9 GeV (See Table. 1). 
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Table 1. Properties and branching ratios of 7r + 7r resonances 



Label % 


Mi 


J p 


m, (MeV) 


Ti (MeV) 


BTZi(Mi -> 7T+7T-) 





M = / (980) 


0+ 


980 


40 ~ 100 


0.52 


0' 


M w = / (1500) 


0+ 


1500 


112 


0.3 


1 


Mi = p(770) 


1- 


770 


151 


1 


V 


M r = p(1700) 


1- 


1700 


240 


0.3 


2 


M 2 = / 2 (1270) 


2+ 


1275 


186 


0.56 


3 


M 3 = D° 


0" 


1865 


1.59 x 10" 9 


1.53 x 10" 3 



physics. 

In Section 2, we present our formalism dealing with B — > ttttIv decays within the 
SM and in its extensions, and the observable asymmetries are considered in Section 3. 
Section 4 contains our numerical results and conclusions. Presented in Appendix are all 
the relevant formulae we use here. 

2 Theoretical Details of Decay Amplitudes 

A. Within the Standard Model 

The decay amplitudes for the processes of Fig. 1, with Mj listed in Table 1, 

B-(p B ) - M.fe, Aj) + W*(q) - n + (p + ) + 7r(p_) + Z"(p,, A,) + P(p„) (1) 
are expressed as 

= -^EE^( r (p^O^)b' Mt |o>(M i ( Pi ,A i )|J iM | J B-(p B )) 

xn j ( SM )(7r + (p + )7r-(p_)||M l (p J ,A l )), (2) 

where Aj = for spin states (/o and £)), A« = ±1, for spin 1 states (p), Aj = ±2, ±1, 
for spin 2 states (/ 2 ), and A; is the lepton helicity, ±|. 
The leptonic current is 

j" = ^y(i-75M, (3) 

and for the hadronic currents we have 

^ = ^(1-75)^6, V t = V ub , Ci = -j= for label i = (0 , 1 (0 , 2; (4) 
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Figure 1: Diagrams for B — > MjW* — > 7t + tt~1~vi decays. 



and 



J3 = 47 M (1 - 75)^, V 3 = K*, c 3 = 1 for label i = 3, (5) 

where q stands for the isospin factor especially due to wu-mesons. We assume that the 
resonance contributions of the intermediate states can be treated by the Breit-Wigner 
form, which is written in the narrow width approximation as 



s M - m\ + iniiTi 



M^m) = I — > ( 6 ) 



where s M = (p+ +P-) 2 and the m/s and r,'s are the masses and widths of the resonances 
respectively (See Table 1). For the decay parts of the resonances we use Q 



<7r + (p + )7r-(p_)||M^,A ? ;)> = 0^^U(0*,n (7) 

where Y"j m (#,</>) are the J = I spherical harmonics listed in Appendix B, and the angles 
6* and 0* are those of the final state ir~ specified in the Mi rest frame (See Fig. 2c). 
The couplings of Mj to tttt are effectively taken into account by the branching fractions, 

BKi(Mi -> 7T + 7T-). 

In order to obtain the full helicity amplitude of the B — ► Tinlv decay, we first consider 
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(a) B rest frame (b) W* rest frame (c) M, rest frame 

Figure 2: The decay B -> M { W* -► (n + n-)(W) viewed from the (a) B~, (b) VT* and (c) 
Mi rest frames. 

the amplitude of B — > MiWi ||, denoted as A4\' 

Ml = -^VicS'iPuXiMpJlf^iMiiPi, *i)\Jip\B-(p B )). (8) 

We express the matrix elements -M^ into the following form: 

Ml = %V lCl Y,rix w L X x l w H^ (9) 
V 2 x w 

where for the decays B — ► MiW* and PU* — ► Zz/, respectively, 

= ^<^i(R.^)iJri5"(Pfl)>, 

= e^<r( w ,A0^b' M |o>, (10) 

in terms of the polarization vectors ew = e(g, A^) of the virtual W . These e^'s satisfy 
the relation 

-<r = E^e^, (ii) 

where the summation is over the helicities Xw — ±1> 0, s of the virtual W, with the metric 
V± = Vo = -ijs = 1- 

We evaluate the leptonic amplitude LJ in the rest frame of the virtual W (See 
Fig. 2b) with the z-axis chosen along the Mj direction, and the x—z plane chosen as the 
virtual W decay plane, with {pi) x > 0. Using the 2-component spinor technique and 



polarization vectors given in Appendix B, we find 



L± = 2\Jq 2 vd±, L = —2^J q 2 vd , L s =0, 

L± = ±2m l vd , Lq = V2miv(d + - L+ = -2m,tJ, (12) 



where 



1 ± cos 0i 



1 d± — p — , and d = sin#j. (13) 



Here we show only the sign of A; as a superscript on L. Note that the L + amplitudes are 
proportional to the lepton mass m h and the scalar amplitude L~ vanishes due to angular 
momentum conservation. 

For B — > Mi transitions through the weak charged current 

J? = Vt - Af, (14) 

the most general forms of matrix elements are 

for /o(0 + ) states : 

(fo(Pi)K\B(p B )) = 0, 

(fo{Pi)\A p \B(p B )) = u +((1 2 )(Pb + Pi)» + U-(q 2 )(p B ~ Pi)p, 
for p(l~) states : 

(p(p i ,e 1 )\V IJ ,\B(p B )) = ig{q 2 )e ilvpa e* l u (p B + pi) p (p B - pi) a , 
(p(p i ,e 1 )\A lx \B(p B )) = f(q 2 )e* lfl + a + (q 2 )(el ■ p B )(p B + Pi)^ 

+a_(g 2 )(e* • p B )(p B - 

for /2(2 + ) states : 

(f2(pi,e 2 )\V tM \Bi PB )) = ih(q 2 )e, uXp e* 2 » a p Ba (p B +p z )\p B -p z y, 

(/ 2 fe,e 2 )|A^5(^)) = A;(g 2 ) e ^ + 6 + (g 2 )( e L /3 ^)(pB+^)M 

+b-(q 2 M*pP a BPB)(PB-Pi)„; 

for D(0~) states : 

(D( P i)\V,\B(p B )) = f + (q 2 )(p B + Pl ), + f-(q 2 )(PB-p l ),, 
(D(p t )\A ll \B(p B ))=0, (15) 
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where e\ and e 2 are the polarization vectors of the spin 1 and spin 2 states, respectively. 
Using the above expressions and the polarization vectors given in Appendix B, we find 
non-zero B — > M{W* amplitudes are 

for i = 0, H° Xw =S° Xw1 

S° = -u + {q 2 )^f, S° = - ( M+ (g 2 ) (m ^ J + u-tf)^) , (16) 



for ^ = l i ' ) , H*=V&, 

1 



V° 



/(g 2 )«-s M -g 2 ) + a + ( g 2 )g + g. 



V? 



f(q 2 )+a + (q 2 )(m 2 B -s M ) + a.(q 2 )q 2 



(17) 



for i = 2, Hj? w =T>«, 

1 



rp±l 



2\/6 








" 2v/2 1 




1 





k{q 2 ){m 2 B -s M -q 2 ) + b + {q 2 )Q + Q. 



for i = 3, H°=P° 



2^6 s M V? 
= P° 



[% 2 )T% 2 )vW?-]' 

/,( ? 2 )+6 + (? 2 )(^-«m)+&-(? 2 )? 2 



where 



Here 



P ° = / + (g 2 )^±^, P s ° = / + (g 2 ) ( ™* J* J + /-(? 2 )aM 



Q± = (rn B ± v 7 ^) 2 - 9 2 - 



(18) 



(19) 



(20) 



(21) 



Q+Q~ = A(m|,s M ,g 2 ) 

gives the triangle function A(a, 6, c) = a 2 + b 2 + c 2 — 2(afe + 6c + ca). Combining all the 
formulae, we can write the full helicity amplitudes of B~ — > ix+ix~l~v decays as 



A=0,s 
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A=±l 



(22) 



where 



n /0 = ^=(y/Bn^u + y/Bn^Uo> 



n D = JBTZ 3 U 3 , (23) 



and 



t = ^r- (24) 

Vub 

Note that we use Y$ for the pseudo scalar meson D, which is actually the same quantity 
as Y® = l/\/47r for the scalar meson f except that it changes sign under the parity 
transformation. Concerning the parametrization of £, other CKM factors, such as V* d 
from D° — > 7T + 7r _ decay, are already included in its branching fraction calculation. And 
because we use implicitly Wolfenstein parametrization J7| for CKM matrix, in which the 
complex phases are approximately in the elements Vm and V u b only, the imaginary part 
of £ here comes only from the element V u b- 

The differential partial width of interest can be expressed as 

2m B ~* 25Qir 3 m%q 2 
where the 4 body phase space <i$ 4 is 



<i$4 = ds M ■ dq 2 ■ dcosQ* ■ dcosOi ■ d(f>* . (26) 
Kinematically allowed regions of the variables are 

Ami < s m < ( m B - m) 2 , 
mf < q 2 < (m B - ^s^) 2 , 
— 1 < cos 9*, cos9i < 1, 

< <p* < 2tt. (27) 

Since the initial B~ system is not CP self-conjugate, any genuine CP-odd observable 
can be constructed only by considering both the B~ decay and its charge-conjugated 
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B + decay, and by identifying the CP relations of their kinematic distributions. Before 
constructing possible CP-odd asymmetries explicitly, we calculate the decay amplitudes 
for the charge-conjugated process B + — > 7r + 7r~/ + z^. For the charge-conjugated B + decays, 
the amplitudes can be written as 

A Xl = -%EE^<i + (P«,A 1 )i/(p,)|jlO><]i? , i (p i) A i )|4 i |B + (pB)> 



i Xi 



V2 



xU t (s M )(7T + (p + )7r-(p_)\\M t ( Pi ,X,)}. 



[28) 



The leptonic amplitudes L$L. are 



L% = -2^/q 2 vd T , L+ = -2^/q 2 vd , Lj 
L± = ±2mivd , Lq = V2rriiv(d + — d_), L~ s 



-2rriiv. 



(29) 



And the transition amplitudes H'S for B + — > MiW* are given by a simple modification 
of the amplitudes H^L of the B~ decays: 



Hx w = Hti9 ^~9,h^ -h, f ± -> -/ ± }. 
Then, we find the full amplitude for B + — > 7r + 7r~/ + z/: 

"C°V0 i t*TT D°V0 



(30) 



where 



'72 



A=0,s 



<V 0) 



A=±l 



(31) 



D A W 


q0 

_ D A W ' 






- _p0 

Aw> 




v 0,s 


- v° F ±x 

— K 0,s> K ±1 


- V Tl 

- V T 1 > 


J 0,s 


— i 0,s' J ±1 





(32) 

It is easy to see that if V u b and V c b are real, the amplitude (|22|) of the _B~ decay and 
II) of the 5 + decay satisfy the CP relation: 



A 



ft) = ricpAT(6*, -0*, ft; F ° -> -y °), 



(33) 



where 9* and <fi* in A Xl are the angles of the final state 7r + , while those in A Xl are for 7r~. 
Then, with a complex weak phase £, dT/d$4 can be decomposed into a CP-even part S 
and a CP-odd part V: 

,/J J(5 + X>). (34) 



2 

The CP-even part S and the CP-odd part T> can be easily identified by making use of the 
CP relation (|33"D between 5~ and B + decay amplitudes, and they are expressed as 

„ rf(r + r) d(r-T) , , 

where T and T are the decay rates for B~ and respectively, and we have used the 
same kinematic variables {s M ,q 2 ,8* ,9i} for the dT/d^^ except for the replacements of 
(f>* -> -<p* and F ° -> -y °, as shown in Eq. (|33|). The CP-even 5 term and the CP-odd 
T> term can be obtained from B T decay probabilities, and their explicit form is listed in 
Appendix C. Note that the CP-odd term is proportional to the imaginary part of the 
parameter £ in Eq. (|24|). 

Before we go further on to the beyond the SM analyses, we note that in addition to the 
resonant tree diagram contributions there are other SM contributions through annihilation 
diagrams and electroweak penguin diagrams, which are relevant for nonresonant case. As 
written in Section 1, we consider only resonant contributions by assuming nonresonant 
contributions can be separated through data analyses. 

B. With complex scalar couplings 

Next we consider CP violation effects in extensions of the SM, where we extend the 
virtual jy-exchange part in Fig. 1 by including an additional scalar interaction with com- 
plex couplings. First we describe the formalism in a model independent way, but later we 
consider specific models such as multi-Higgs doublet models and scalar-leptoquark mod- 
els. In this case CP-violating phases can be generated through the interference between 
jy-exchange diagrams and scalar exchange diagrams with complex couplings. 

The decay amplitudes for B — > 7T + 7T~l~i?i are expressed as 

A x ' = -K fe %^EEf( r b^^(^)ir t |0)(M l (p,,,A i )|J M | J B-( Pi? )) 
V2V2 i K L 
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H(r(pi,h)Hp»Mu)(Mi(Pi,\i)\J s \B-(p B )) 

xn j ( SM )(7r + (p + )7r-(p_)||M i fe,A,)), (30) 



where 



r = ^(1 - 75)^ = V» - A", (37) 

and their corresponding scalar currents are 

j a = V>„(1 - 75)^, ^ = - 75)^6, (38) 

the additional factor l/v2 comes from the isospin factor as mentioned earlier. Here the 
parameter (, which parameterizes contributions from physics beyond the SM, is in general 
a complex number. And as explained earlier, in order to exclude any possible CP violation 
effects induced within the SM, we only include the lowest three states p(770), / (980) and 
/ 2 (1270) as intermediate states. By using the Dirac equation for the leptonic current, 
%j tl = mrfs, the amplitude can be written as 

A x ' = -K b ^^^(rfe,AO^)|^ t |0)(M i ( Pi ,A i )|^| J B-(p B )) 

xn j ( SM )(7r + (p + )7r-(p_)||M l (p l ,A,)), (39) 

where the effective hadronic current is defined as 

Q = J +(^LJ S . (40) 
mi 

In this case the amplitudes of B — > M^lv have the same form as the previous SM 
case (^) except for the modification in the hadronic current part due to the additional 
scalar current: 

Mi = ^V ub J2 VXw L% w Ht, (41) 
where Ji^ stands for the hadronic amplitudes modified by the scalar current J s . Using 



the equation of motion for u and b quarks, we get within the on-shell approximation 



Js = (Pb - Pu) 



V, A, 
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(42) 



Later for numerical calculations, we use the approximation, (p^ — p^) « (p^ — p^.) = 
g M , which is assumed in quark model calculations of form factors ||17|| . After explicit 
calculation, we find that the additional scalar current modifies only the scalar component 
ofH x x - i-e- 

K = (1-C)H° S , 
and Ht = Ht for X w = 0, ± 1, (43) 



where 

C = ( q \ X . (44) 
mi(m b + m u ) 

In this case, dT /g?$4 also can be decomposed into a CP-even part S and a CP-odd part 
V: 

1,1 Vs + V). (45) 



d$ 4 2 

Their explicit form is listed in Appendix C. Note that the CP-odd term is proportional 
to the imaginary part of the parameter ( and the lepton mass m^. Therefore, we have to 
consider massive leptonic (fi or r) decays. 

As specific extensions of the SM, we consider four types of scalar-exchange models 
which preserve the symmetries of the SM H: One of them is the multi-Higgs-doublet 



(MHD) model |§ and the other three models are scalar-leptoquark (SLQ) models |T0 



]T|. The authors of Ref. [12| investigated CP violations in r decay processes with these 
extended models. We follow their description and make it to be appropriate for our 
analysis. 

In the MHD model CP violation can arise in the charged Higgs sector with more than 
two Higgs doublets [ff3l and when not all the charged scalars are degenerate. As in most 
previous phenomenological analyses, we assume that all but the lightest of the charged 
scalars effectively decouple from fermions. The effective Lagrangian of the MHD model 
contributing to the decay B — > ttttWi is then given at energies considerably low compared 
to M H by 

C MHD = 2V2G F V ub ^\m b XZ*(u L b R ) + m u YZ*{u R b L )} (l R u L ), (46) 
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where X, Y and Z are complex coupling constants which can be expressed in terms of 
the charged Higgs mixing matrix elements. From the effective Lagrangian, we obtain for 
the MHD CP-violation parameter lm(( MHD ), 

MC MHD ) = ^{MXZ*) - (^)Im(yZ*)}. (47) 

The constraints on the CP-violation parameter ( ]4"T| ) depend upon the values chosen for 
the u and b quark masses. In the present work, we use (See Appendix A) 

m u = 0.33 GeV, m b = 5.12 GeV. (48) 

In the MHD model the strongest constraint on lm(XZ*) comes from the measurement 
of the branching ratio B(b — > Xru T ), which actually gives a constraint on \XZ\. For 
M H < 440 GeV, the bound on lm(XZ*) is given by 

lm(XZ*) < \XZ\ < 0.23M^ GeV" 2 . (49) 

On the other hand, the bound on Im(YZ*) is mainly given by K + — > 7i + uu. The present 
bound || is 

lm(YZ*) < \YZ\ < 110. (50) 
Combining the above bounds, we obtain the following bounds on lm(( MHD ) as 

|Im(C MH J| < 2.06 for r family, 

|Im(C MH J| < 0.12 for // family. (51) 



On the other hand, the effective Lagrangians for the three SLQ models M, 10 



con- 



tributing to the decay B — > ttttIu are written in the form, after a few Fierz rearrangements: 

+ /i.e., 



c 1 

SLQ 



SLQ 





















(5 


2M? 

02 




Z 3j Z tj 


(5 


2 < 





(b L u R )(v lL l R ) + ^(&lct M! J u R )(vi L o ^Jjt) 



(hu R )(l R v? L ) + \{b L a^u R W R a^v c lL ) 
(b L ^u L )(l c L ^vf L )+h.c., 

^ 1V1 4>2 

Cg = - Z ^(h^u L ){llY^ L ) + h.c. , (52) 
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where j — 2, 3 for I = /j,t, respectively and the coupling constants a£- , j/^ and % are in 
general complex so that CP is violated in the scalar-fermion Yukawa interaction terms. 
The superscript c in the Lagrangians £'* LQ and C* denotes charge conjugation, i.e. 
ip c RL = ^7°7 2 V^l in the chiral representation. Then we find that the size of the SLQ 
model CP-violation effects is dictated by the CP-odd parameters 



Im(C 



i N _ lm[x 3j x'^} 



W 4V2G F V ub Ml 



,/* i 



mO=°- ( 53 ) 

Although there are at present no direct constraints on the SLQ model CP-odd parameters 
in (|53"|), a rough constraint to the parameters can be provided by the assumption |14| that 



\ x ij 



\xij\ and \y[j\ ~ \yij\, that is to say, the leptoquark couplings to quarks and 
leptons belonging to the same generation are of a similar size; then the experimental 
upper bounds from BB mixing for r family, B — > fip,X decay for \x in Model I, and 



B — > IvX for t together with the V u b measurement for \i in Model II yield |14 



|Im(4 Q )| < 2.76, \hn{C Q )\ < 18.4 for r family, 
MOl < 0-37, |Im(C" iQ )| < 1.84 for n family. (54) 
Based on the constraints (O) and fl5"4]) to the CP-odd parameters, we quantitatively 



estimate the number of B — > 7r + 7r / ^ decays to detect CP violation for the maximally- 
allowed values of the CP-odd parameters. 

3 Observable CP Asymmetries 

An easily constructed CP-odd asymmetry is the rate asymmetry 

r - r 

A = =, (55 

r + r' v ; 

which has been used as a probe of CP violation in Higgs and top quark sectors [ 15] . Here 



T and T are the decay rates for B and B + , respectively. The statistical significance of 
the asymmetry can be computed as 

N--N+ N„-N + 

NsD = 7^=m = 7F=^' (56) 
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where N$d is the number of standard deviations, N± is the number of events predicted in 
B14 decay for B ± meson, N is the number of S-mesons produced, and Br is the branching 
fraction of the relevant B decay mode. For a realistic detection efficiency e, we have to 
rescale the number of events by this parameter, N_+N + — > e(iV„ + iV + ). Taking N$d — 1) 
we obtain the number Nb of the B mesons needed to observe CP violation at 1-er level: 

- UTTF (57) 

Next, we consider the so-called optimal observable. An appropriate real weight func- 
tion w(s M ,q 2 ;9*,9i,<p*) is usually employed to separate the CP-odd V contribution and 
to enhance its analysis power for the CP-odd parameter through the CP-odd quantity: 



(wV) = J [wV] d$ 4 , (5S 



and the analysis power is determined by the parameter, 

(wV) 



(59) 



For the analysis power e, the number Nb of the I?-mesons needed to observe CP violation 
at l-o" level is 

n b = IT— 2 ■ ( 6 °) 
Br ■ e z 

Certainly, it is desirable to find the optimal weight function with the largest analysis 
power. It is known |16] that when the CP-odd contribution to the total rate is relatively 



small, the optimal weight function is approximately given as 



V 

w opt (s M ,s L ;9,9 h (p) = — =>• e opt 



S opt \ (S) ' 

We adopt this optimal weight function in the following numerical analyses 

4 Numerical Results and Conclusions 



^ (61) 



Now we show our numerical results. We use the so-called ISGW predictions [17| for all 
the form factors in B — > Mi transition amplitudes of Eq. (|I5|) . One can find in Ref. Jl7| 
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Table 2. The CP- violating rate asymmetry A and the optimal asymmetry e opt , determined 
within the SM, and the number of charged B meson pairs, N B , needed for detection at 
la level, at reference value Im(£) = 12.5. 

B — > 7T + 7T~li?l 



Mode 


I = 


e 


I = 




I = 


r 


Asym. 


Size(%) 


N B 


Size(%) 




Size(%) 




A 


0.94 x 10" b 


1.37 x 10 iS 


1.71 x 10 -6 


4.16 x 10 i7 


1.14 x 10" e 


1.46 x 10 iS 


^opt 


1.45 x 10~ 2 


5.75 x 10 9 


1.44 x 10" 2 


5.79 x 10 9 


1.11 x 1(T 2 


1.56 x 10 10 



the detailed description of the general formalism and relevant form factors for B — > Xeu e 
after neglecting lepton masses. In Appendix A, we give explicit expressions of form factors 
needed for semileptonic decays with non-zero lepton masses. 

We first consider the case within the SM. Total branching ratio of the B~ — ► M — ► 
7r + 7r _ )ez/ e , Mi = p,fo t 2,D is about 0.8%. It depends on the chosen value for \V u b\, and 
here we adopt the result by CLEO |T8] 



w, 



ub\ 



3.3 ± 0.4 ± 0.7 x 10" 



(62) 



In Table 2, we show the results of B — > Tmlu decays for the two CP-violating asymmetries; 
the rate asymmetry A and the optimal asymmetry e opt . We estimated the number of B- 
meson pairs, N B , needed for detection at la level for maximally- allowed values of CP-odd 
parameters Im(£) in Eq. (|24]). We use the current experimental bound |19| 

V ub 



V, 



cb 



0.08 ±0.02, 



(63) 



which means 



|f | = 12.5 ±3.13. 



(64) 



The results in Table 2 are for the maximal case with Im(f) = |f | = 12.5. Due to the 
large cancelations in the simple rate asymmetry [|J when we integrated over the phase 
space, the optimal observable gives much better result. For example, using the optimal 
observable, we need ~ 10 9 5-meson pairs to detect the maximal CP-odd effect in electron 
mode. CP violation effects in B — ► 7t + tt~Wi decays within the SM are not likely to be 
detected, with O(10 s ) 5-meson pairs to be produced at the asymmetric B factories. One 
may rely on hadronic S-factories of BTeV and LHC-B. 
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Table 3. The CP- violating rate asymmetry A and the optimal asymmetry e opt , determined 
in the extended models, and the number of charged B meson pairs, N B , needed for 
detection at la level, at reference values (a) lm(( MHD ) = 2.06, lm(( SLQ ) = 2.76 and 

Im(C" iQ ) = 18.4 for the B tA decays and (b) Im(C Am J = 0.12, hn(C' SLQ ) = 0.37 and 

Im(£ ) = 1.84 for the B^ decays. 



(a) B — > 7r + 7r tv t mode 



Model 


MHD 


SLQ I 


SLQ II 


Asym. 


Size(%) N B 


Size(%) N B 


Size(%) N B 


A 

^opt 


1.47 x 10" 3 7.63 x 10 11 
16.2 6.23 x 10 3 


2.67 x 10" 3 1.99 x 10 11 
18.2 4.27 x 10 3 


3.62 x 10" 3 7.39 x 10 9 
9.67 1.04 x 10 3 




(b) B — > 7r + 7r fiv^ mode 


Model 


MHD 


SLQ I 


SLQ II 


Asym. 


Size(%) N B 


Size(%) N B 


Size(%) N B 


A 

^opt 


2.61 x 10" 5 1.93 x 10 15 
0.18 3.89 x 10 7 


0.90 x 10" 4 1.59 x 10 14 
0.50 5.13 x 10 6 


3.43 x 10" 4 8.89 x 10 12 
1.48 4.76 x 10 5 



Next we consider the extended model case. In this case, CP violation effects are 
proportional to the lepton mass, and we consider only massive lepton (/i or r) cases. In 
Table 3, we show the results of _B r4 and B^ decays. Here in order to distinguish new 
physics effect from the SM one, we use a cutoff for the invariant mass of the final state 

7T + 7T~ aS 

v ^7 < 1.4 GeV. (65) 

We consider only the lowest three uu states, p(770), / (980) and / 2 (1275) as intermediate 
resonances in Table 1 so that the effects of D meson can not enter, and we can thus 
ensure that the result is solely from new physics. Similarly as in the SM case, we estimate 
the number of B- meson pairs, N B , needed for detection at la level for maximally- allowed 
values of CP-odd parameters Im(£) of Eq. ( [5i"l) and (|54|). We again find the optimal 
observable gives much better results than the simple rate asymmetry. 

The results in Table 3 show that CP violation effects from new physics are readily 
observed in the forthcoming asymmetric .B-factories, by using optimal observables. As 
expected, B T ^ decay modes give better results than B^ cases for the MHD model, where 
the CP-odd parameter itself is proportional to the lepton mass. For example, the current 



16 



bounds in the MHD model 



Im(C MH J = 2.06 (0.12) for r (/i) 



directly result from the lepton mass dependence. But there is no such dependence in the 
SLQ models. The current numerical values of CP-odd parameters in the SLQ models, 



are just from different experimental bounds. Therefore, the smaller CP-odd value for p 
family is a consequence of the fact that the current experimental constraints on the muon 
mode are more stringent. And B r4 decay modes would provide more stringent constraints 
to all the extended models that we have considered. 

In conclusion, we have investigated direct CP violations from physics beyond the SM as 
well as within the SM through semileptonic decays: B ± — > ix+ix^l^vi. Within the SM, 
CP violation could be generated through interference between resonances with different 
quark flavors, that is, with different CKM matrix elements. We included uu state mesons 
(p, fo and / 2 ) and D meson as intermediate resonances which decay to it + -n~. Using 
optimal observables, we found (9(10 9 ) _B-meson pairs are needed to probe CP violation 
effects at la level for the current maximal value of Im(£) = IKa/V^I = 12.5. We have 
also investigated CP violation effects in extensions of the SM. We considered multi-Higgs 
doublet model and scalar-leptoquark models. Here CP violation is implemented through 
interference between VU-exchange diagrams and scalar-exchange diagrams with complex 
couplings in the extended models. We calculated the CP-odd rate asymmetry and the 
optimal asymmetry for £? r4 and B^ decay modes. We found that the optimal asymmetries 
for both modes are sizable and can be detected at la level with about 10 3 -10 7 5-meson 
pairs, for maximally-allowed values of CP-odd parameters. Since ~10 8 5-meson pairs are 
expected to be produced yearly at the asymmetric B factories, one could easily investigate 
CP-violation effects in these decay modes to extract much more stringent constraints 
on CP-odd parameters, lm(( MHD ) and Im(<^"). 




2.76 (0.37) 

18.4 (1.84) for r (p), 
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Appendix 



A Form factors 



Form factors in Eq. (fLq) within ISGW model [17j are 

m u m b m q 



u + (q 2 ) = -F 5 (q*;f ) 



\/6/3 B m/ /i- 



/ 2\ t-i / 2 /• ^ u (m B + m fo ) 
u+[q ) = F 5 (q ;/ )- 



y/6(3 B m fo 



(66) 



a+(g 2 ) 
a_(g 2 ) 



2^ ^s(g 2 ;p) 



m g 2[i-m p f3 Bp 



f(q 2 ) = 2m B F 3 (q 2 ;p), 



F 3 (q 2 ;p) 
2m n 



1 + 



m u fP 2 B-P 2 p \ rn 2 J 2 p 



m b \f3% + (3 2 p ) 4p_m B f3 Bp 



2^ ^(<? 2 ;p) 



2m , 



m b \ m q f3 Bp 



m 2 u B 



2 pa 
p 



4p + m B 3 Bp 



(67) 



Kq 2 ) = F 5 (q 2 ;f 2 ) 
b + {q 2 ) 
b-{q 2 ) 



m. 



2V2m B (3 B 
F 5 {q 2 ] h)m u 



1 

m n 



m u l3 B 



2V2f3 B m h 



m b 



2//_m /2 /% 2 



k(q 2 ) 



V2^F 5 (q 2 ;f 2 )., 



m u (3 2 fo mlm b (3f 2 



•j, _ F 5 (q 2 ;f 2 )m u 



2\[2fi B m h m b 



1 + 



rn 2 u m b {3% 



m q m B (3 2 Bh Afx + m B P Bf . 



(6* 



where 



And 



U(q 2 ) = F 3 (q 2 ;D) 
f.(q 2 )=F 3 (q 2 ;D) 



1 _|_ _ m hm q m u (5 2 B 
2/x_ 4p + p^m D (3 BD 
/ 1 

1 - (me + rh D ) 



m u (3 



2m q 4p: + m D (3 BD J 



?bx ~( >' /'±- ( — ± — 



1 1 
— ± — 

m q m b 



-i 



^n(g 2 ;X) 



/%\ 1/2 (Wx 
\m B ) 1/3; 



n/2 



exp 



/ m 2 \ g m - g 2 
\4m B m x ) K 2 f3 2 BX 



(69) 



(70) 



(71) 



where relativistic compensation factor k — 0.7, and q m is the maximum value of q 



2. 



q m = (rn B - y/s 
19 



2 

Ml 1 



(72) 



and m q is m u for uu state mesons and m c for .D-mesons. The numerical values of fix m 
GeV unit are 

Pb = 0.41, (3 D = 0.39, p fo = 0.27, f3 p = 0.31, (3 h = 0.27, (73) 
and quark masses in GeV unit are 

m u = 0.33, m c = 1.82, m b = 5.12. (74) 
The so-called mock meson masses rhx are defined as 

rh B =m b + m u , m D =m c + m u , m p j ,f 2 = 2m u . (75) 

B Kinematics 

Spherical harmonics 



Y° 



Y° 

I 2 



47r' 




47T 

5~ 3 
4tT2 



cos^, Ff 



±i 




sinfle^, 



7z(^cos 2 ^- h, F 2 ±1 = =Ft/ ^ sin^ cos fle* 
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8tt 



Polarization vectors 



(76) 



In the I? rest frame, where the coordinates are chosen such that the z-axis is along the 
Mi momentum and the charged lepton momentum is in the x-z plane with positive x- 
component (cf. Fig. 2a), the polarization vectors for the virtual W are 

1 



e(q,±r 
e(q,0r 
e(q,sr 



T — (0,1,^,0), 
-±=(p M ,0,0,-q°), 



(77) 
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and the polarization states of the spin 1 mesons are 

e(±l)'* = =F-^(0,l,±i,0), e(py = -^=(p M ,0,0,E M ), (7* 



where p M = \ / Q+Q-/2mB with Q± defined in Eq. (|20|) , and E M = (m% + s M — q 2 )/2rriB- 
And for the spin 2 meson we get 

e{±2Y v = e^(±l)e^(±l), 

e(±ir = ^[^(±1K(0)+^(0K(±1)], 

e(0r = ^ [e p (+l)£"(-l) + ^(-1K(+1)] + ^(0K(0). (79) 
In the W rest frame the polarization states of the virtual W are 

e(g,±) M = T^(0,1,T*,0), 
e( 9 ,0)" = (0,0,0,-1), 

e(g,s)" = ~^q» = (1,0,0,0). (80) 

C CP-even and CP-odd quantities 

Within the SM 

The CP-even quantity S is 

S = 2C(q 2 ,s M )E, (81) 

with 

s = (^o ^o ) 2 !^,! 2 + (^o n ^o°) 2 iei 2 |nz,| 2 + |{v-)n p | 2 + |{T->n /2 | 2 
+2(L -^r Q )Re(n /0 n;(y-)* + n /0 n /2 (T-)*) + 2Re(n p n} 2 (y->(T-)*) 
+2(L -p r )Re(0[(£ 5 °y °)Re(n D n} ) + Re{u D w p (v-y + n D n} 2 (T-)*)] 

+(L+SX° - LtS%°m \ 2 + (L+P°Y ° ~ LtPXWW 

+ l n p| 2 K^ + ) - L t V s Y l\ 2 + l n / 2 | 2 K T+ ) - L t T s Y 2\ 2 

+2(l+5 ° - L+5 s )r [-(L+K F 1 )Re(n /0 n;) + Re(u fo n;(v + y) 
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-(L+T s °F 2 °)Re(n /0 n* 2 ) + Re(U fo U%(T + )*)} 

+2Re(n p n* 2 (\/ + )(T + )*) + 2(LiV s %°)(LjT°Y 2 °)Mn p n} 2 ) 

-2(L+T s °r 2 °)Re(n p n} 2 (\/ + )) - 2(L+\/°F 1 °)Re(n p n} 2 (T + )*) 

+2(L+P° - LtP?)Y °RemLtS° - S° s )Y°Re(n D n%) - (L+ W)Re(n D n;) 

+Re(U D U* p (V + )*) - (L+T s °F 2 °)Re(n D n} 2 ) + Re(U D U* h (T + )*)}, (82) 

and the CP-odd quantity V is 

V = 2Im(0C(q\s M )A, (83) 

with 

A = -2(L -P %°)[(L -SM^ 

-2(z+p ° - LtP°)Y °{(L+s° - Lfs° s )Y °im(n D n%) - (L+K°n°)im(n D n;) 
+im(n D n p (\/ + >*) - (L+Tfr 2 °)im(n D n} 2 ) + im(u D u} 2 (T + y)\, (84) 



where 



(V±) = £ LtVfYf, (T ± ) = £ LtnY 2 \ (85) 

i=0,±l i=0,±l 

and the overall function C(q 2 , s M ) is given by 

r(n 2 s ,2^1 i (g 2 -mf)yg + g_ 

With a complex scalar coupling 

The CP-even quantity S is 

S = 2C{q\s M )^ (87) 

with 

s = (Vo%°) 2 |n/ol 2 + ir)n/ + |(r)n /2 | 2 

+2(L -5 r °)Re(n /0 n;(\/-)* + n /o rr} 2 <r->*) + 2Re(n p rr /2 <y-><r->*) 

+|n /0 | 2 |L ^ % -(i-c')^ s %T 
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+ |n p | 2 [|(^ + )| 2 + (L s + W) 2 |l - CI 2 - 2(LfV s °Y?)Re((V+))Re(l - C')] 
+ |n /2 | 2 [|(T + )| 2 + (L+T s °y 2 °) 2 |l - Cf - 2(L+T°Y 2 °)Re((T+))Re(l - C')] 
+2Re(n /o n;)[(L +S - L+ S° s )Y °Re((V + )) - W)Re(l - C') 

+ (L+SX°)Re((^ + ))Re(C') + (^ s °y o °)(L+ W)|l - C'| 2 ] 
+2Im(n /o n;)Im((V+>)[(L+S ° - L+S?)ltf + (L+^y o )Re(C')] 
+2Re(n /0 IT /2 )[(L +S - LiS° s )Y °Re((T+}) - (^«)(L+T°y 2 °)Re(l - 
+ (^OT)Re((T + ))Re(C') + (^5 s °r o )(L+T s r 2 °)|l - Cf] 
+2Im(n /0 n* 2 )Im((T + ))[(L + S ° - L+S?)y o ° + (L+^)He(0] 
+2Re(n p n* 2 )[Re((V + )(T + )*) - (L+T°F 2 )Re((V + )) + (L+TjY 2 °)Re(<V + ))Re(0 
-(L+K°n°)Re((T + ))Re(l - C') + - C'| 2 ] 

-2Im(n p n} 2 )[Im((V + )(T + )*) - (L+T°r 2 °)Im((y + )) + (L+T°y 2 °)Im((V+))Re(C') 
+ (L+K^Im((T+>)Re(l -(')], 



and the CP-odd quantity V is 

P = 2Im(C / )C(g 2 , SM )A, (89) 

with 

a = 2[im((\/ + )){(L+K°F 1 °)|n p | 2 + (L+s s °r )Re(n /o np + (^T°r 2 °)Re(n p n} 2 )} 
+im((T + )){(L+r°r 2 °)|n /2 | 2 + (Lt s°X)Re(u fo u%) + (L+K°n°)Re(n p n} 2 )} 
+Re((v + )){(L+T°y 2 °)im(n p n} 2 ) - (LXX°)Mn /o n;)} 
-Re((T+)){(L+v; r 1 )im(n p n} 2 ) + (L+s^)im(n*m, )} 



+(^^o n )(^K n )im(n /o n;) + (L + 5 F °)(L+T°y 2 °)im(n /o n 



(90) 



Note that since every term in A of Eq. (90) contains square terms of Lf which are 
proportional to mi (see Eq. (0)), the CP-odd quantity T> of Eq. (89) is proportional to 
lepton mass due to the definition of (' (see Eq. (fH|)). 
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